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Abstract 

In this article, we give a brief informal introduction to Malliavin Calculus for newcomers. 
We apply these ideas to the simulation of Greeks in Finance. First to European-type options 
where formulas can be computed explicitly and therefore can serve as testing ground. Later 
we study the case of Asian options where close formulas are not available. The Greeks are 
computed through Monte Carlo simulation. 

1 Introduction 

Malliavin Calculus is an area of research which for many years has been considered highly 
theoretical and technical from the mathematical point of view. In recent years it has become 
clear that there are various applications of Malliavin Calculus as far as the integration by parts 
formula is concerned. Nevertheless it is still considered by the general mathematical audience 
and practitioners a field where is hard to grasp the basic ideas or to obtain new contributions. 

In this article we present an informal introduction to Malliavin Calculus which we hope can 
open the area to practitioners. For serious mathematical approaches to the topic we refer the 
readers to the authoritative books on the matter, e.g. [||, ||, We have tried to take the spirit 
of the issues to try to explain in simple terms the elements of the theory. 

2 Malliavin calculus 

The most common concise way Malliavin Calculus is presented in a research paper is as follows. 
Let W = {VFj}tg[o,i] be a standard one-dimensional Brownian motion defined on a complete 
probability space {Q,J^,P). Assume = {J^t}te[o,i] is generated by W. Let S be the space 
of random variables of the form F = /(Wj^, . . . , Wt„), where / is smooth. For F € S, DfF = 
Er=i m-fi^ti^- ■ ■ ' WtJl[o,u]it). For k eZ+,p>l, let D*='P be the completion of S with the 
respect to the norm 

||F||,,, = {E[\Fn + /'.../' \Di^,...,s,F\'ds^ . . . ds,r/'])'/^, 

~[Jo Jo 

where D{^^ ^.F = Dt^...Dt^F. We let ||F||o,p = {E[FP]y/P = \\F\\p and = nfc,pD*^'P. For 
processes u = {?it}tg[o.i] on {Q,J^,P), B^fj-jQ -^^^ is defined as B'^'P but with norm ||^i||fc,p,L2([o,i]) = 

(^[ii^iii^do,!])] + ^[(E •=! /; . . . /; \\Dk,... ,sAhm])dsi . . . ds,r/'])'/p. 

We denote by D*(u) the Skorokhod integral or the adjoint operator of D. This adjoint 
operator behaves like a stochastic integral. In fact, if ut is J^t adapted, then D*{u) = utdWt, 



the Ito integral of u; see e.g. p. Here we write D*{u) = /J utdWt, even if ut is not J^t adapted. 
There are other anticipating integrals that have some relationship with this one as e.g. the 
Ogawa symmetric integral. Of the formulas we will use, the following are worth mentioning, 



FutdWt = F / utdWi 



{DtF)utdt, 



for F G ©1-2 and ^[F^ u'fdt] < oo — see e.g. ||— ; and 



E 



{DtF)utdt 



E[FD*{u)]. 



As a byproduct of all the above formulas one obtains the integration by parts formula. For 
this, we say that F is smooth if F G D°°. For a real random variable F G B^'-^, we denote by V'F 
the Malliavin covariance matrix associated with F. That is, V'F =< DF,DF >l2[o,i]xR- One 
says that the random variable is non-degenerate if F G B°° and the matrix ipp is invertible a.s. 
and (detV'F)""^ G r\p>iU'{n). 

The integration by parts formula of Malliavin Calculus can be briefly described as follows. 
Suppose that F is a non-degenerate random variable and G G D°°. Then for any function 
g G C^(M'') and a finite sequence of multi-indexes /3, we have that there exists a random 
variable H^{F, G) so that 



E[g^iF)G] = E[g{F)H^{F,G)] with 



//^(F, G) < C(n,p,/3) \\det{^Pp)-%, \\F\\l, \\G\\,,^,, , 

n,p ^ 

for some constants C{n,p,(3), a, b, d, p', a', b', d' and (3 G U„>i{l, g}". Here g^ denotes the 
high order derivative of order l{f3), the length of the multi-index and whose partial derivatives 
are taken according the index vector (3. 

A gentler introduction may say that the idea behind the operator D is to differentiate a 
random variable with respect to the underlying noise being this generated by the Wiener process 
W . Therefore heuristically one may think that D^" = " q(^^\y ) ' ^^^^ ^^^^ ™ mind one can guess 
how to differentiate various random variables. Some examples are 

DtWt = 1, 
Dtf{Wt) = nWt), and 

Dsl^j^ f{Wu)dW^ = nW^)dWu + f{Ws). 

Here / is a function. A way to understand any integration by parts formula is through the 
following general definition. 

Definition We will say that given two random variables X and Y , the integration by parts 
formula is valid if for any smooth function / with bounded derivatives we have that 

E[nX)Y]=E[f{X)Hl 

for some random variable H = H{X, Y). 

One can deduce an integration by parts formula through the duality principle @. That is, 
let Z = f{X). Then using the chain rule we have 

D,Z = f'{X)DsX. 



2 



From here we multiply the above by Yh{s) where his a process to be chosen appropiately. Then 



DsZYh{s) = f'{X)DsXYh{s). 
Integrating this for s G [0, 1], we have that 

[ DsZYh{s)ds = [ f'{X)DsXYh{s)ds = f'iX)Y [ h{s)DsXds, then 
Jo Jo Jo 



^ YDsZh{s) 
Jq h{v)DyXdv 



ds = f'iX)Y, therefore 



with 



E <DZ,u>L2^o,i]=E[f'{X)Y] , 
Yh{s) 



jQh{v)D^Xdv' 



Finally, we have that if D* is the adjoint operator of D — see equation (|2|) — , then 

E[ZD*{u)] = E [f'{X)Y] , and 



E 



f{X)D* 



Yh{.) 



J^h{v)D.,Xdv ^ 
This also means that in particular for /i = 1 we have that 



E [f{X)Y] 



Y 



H = H(X, Y) = D* \ 

\j^D,Xdv^ 



(3) 



If one has higher order derivatives then one has to repeat this procedure iteratively. The use 
of the norms in the spaces D"'^ is necessary in order to prove that the above expectations are 
finite (in particular the ones related to H). Note that the integral h{v)DvXdv should not be 
degenerate with probability one. Otherwise the above argument is bound to fail. The process h 
that appears in this calculation is a parameter process that can be chosen so as to obtain this 
non-degeneracy. In the particular case that h{v) = D^X one obtains the so called Malliavin 
covariance matrix. 

In conclusion one can build different integrations by parts formulas depending on how we 
choose this process h. In the next section we use this formula in order to apply it to a concrete 
problem in Finance. 



3 Greeks in Finance 

European options are contracts that are signed between two parties (usually a bank and a 
customer) that allows to obtain certain monetary benefits if the price of certain asset fall above 
(call option) or below (put option) a certain fixed value, the strike price, at a certain fixed date, 
the expiration time. A Greek is the derivative of an option price with respect to a parameter. In 
general, let X = X{a) be a random variable that depends on a parameter a. Suppose that the 
option price is computed through a payoff function in the following form V{a) = E[^{X{a)^ a)] 
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where $ is generally non-smooth. A Greek is therefore a measure of the sensibility of this price 
with respect to its parameters. In particular, it could serve to prevent future dangers in the 
position of a company holding these options. The problem of computing Greeks in Finance has 
been studied by various authors: @, |, |, 0, |, among others. Let us take a clooser look at the 
problem. If the Leibnitz rule of interchange between expectation and differentiation were true 
then we would have 



dV{a) _ dE[^{X{a),a)] _^ 
da da 



da da 



(4) 



When the above expression does not have a close formula then one may start thinking in per- 
forming some Monte Carlo simulations in order to approximate the above quantity. If $ is 
somewhat regular then we can use the last expression above to do this: we shall call this proce- 
dure the "direct method" . Unfortunately in various cases $ is not differentiable. Then one can 
resort to the use of the middle expression above to generate what is known as finite difference 
method. This method has been somewhat successful in the recent past but we would like to 
discuss here the application of the integration by parts pesented before in order to compute 
these derivatives. 



4 The European-style options 

We shall illustrate how this procedure works by choosing a very special subset of the large family 
of Financial derivatives: what we have called European- style options. In this class of derivatives 
we will find all the options whose payoff function depends only on the value of the underlying 
at the expiration time T, which is previously fixed. Examples of the European-type options 
are the vanillas — the more classical European calls and puts — , or the binaries — the so called 
"cash-or-nothing" options — , among others. These options will differ, for instance, from the 
American- style options, where the execution time is not fixed but belongs to an interval; and 
also from the Asian- style options where the payoff depends on some average of the value of the 
asset in a given period of time. We will return on this topic afterwards. 

The interest of the European-style options is that they are a class of derivatives whose Greeks 
can be computed in closed form for particular classes of payoff functions. The reason, as we will 
show, is that we explicitly know the probability density function of the random variable involved, 
St, whereas in other scenarios this is not true. This peculiarity provides us with a framework 
where we can easily test how Malliavin Calculus applies to the computation of Greeks. Later, 
we will also make a comment on a case where this closed formulas are not available and where 
this technique may prove useful. 

4.1 The Malliavin expressions 

Let us start deriving the formal expressions for some of the Greeks we shall deal with. 

First we assume that our underlying asset S is described by a geometric Brownian motion 
under the risk neutral probability P: 

St = So+ f rSsds + f aSs dWs, (5) 
Jo Jo 

where r is the interest rate and a is the volatility. This model is one of the models typically 
used to describe stock prices or stock indices. 

Second, from the previous arguments it follows that X{a) must be in general a functional of 
S. In the case of European-type options, X{a) = St and from @: 

5r = 5oe^^^+"^^>, (6) 
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where {Wt}t(z[o^T] is the Wiener process, and is just r — cr^/2. Expression (|6|) is involved in aU 
the following derivations. 

Now we can compute Delta, A, the first partial derivative of the (discounted) expected 
outcome of the option, with respect to the present value of the asset: 



A 



_d_ 

dSr 



-rT 



Sn 



-E 



dSr 



-rT 



5n 



-E [$'(St)5t] . 



Now we may perform the integration by parts applying the formula given in @, 



-rT 



A 



So 



-E 



^{St)D* 



S' 



T 



(7) 



which removes the derivative of ^ from the expectation. 

The integral term appearing in the last expression will appear many times along our expo- 
sition. In order to compute it we must remember the rules of the stochastic derivative stated 
above: 

DuSt = uStDuWt = aST^uKT, 

and then 

DuSt = tTSt. (8) 
Then we are able to perform the stochastic integral in (|^, 



D* 



St 



D* 



St 



D* 



1 

of 



Wt 



y/p DySTdv ^ 

with the help of equation applied to -F = Then the expression for A reads. 



A = E 



SqctT 



(9) 



Let us move now into a new Greek: Vega. It measures how sensitive is the price of the option 
when the volatility changes, 



V = ^E[e-^^<^{ST)] = e-'^E 
oa 



e-''^E[ST{WT - aT)<^'{ST)]. 



We invoke again the recipe in Section 2 and thus we can withdraw the derivative form ^{St), 



V = e-'-^E 



<^>{St)D 



, I St{Wt-(jT) 



\ Jq DySTdv 

where we have used the expression (M). So the computation we must face is 



e-'-^E 



HSt)D* - 1 



D* (^^-l^=^D*iWT)-WT. 



Here a new instance of stochastic integral appears, D*{Wt)- The rule which we must take into 
account in order to solve the problem is again in (|^), with F = Wt, 



D*{Wj 



DsWTds = W^-T, 
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what lead us finally to this expression for V, 



V = E 



-rT 



wl 



Wt 



a 



$(5t) 



(10) 



The last example we will present here is one involving a second derivative: Gamma. T inform us 
on the second order dependence of the price of the option on the actual value of the underlying, 



-rT 



q2 



After a first integration by parts we obtain, 



o2 
•^0 



-E 



The stochastic integral may be simplified using once more formula (|l]) on F = leading to 



o2 



-E 



^ aT aT ^ ^ aT 







^ (Wt 
DsSTds = StI^-1 



Afterwards we can perform the second integration by parts, yielding: 



-rT 



T 



o2 



-E 



^'{St)St!^^-1 



-rT 



o2 



-E 



'^{St)D* 



S' 



T 



J^D.StcIv I 



The stochastic integral is now slightly cumbersome, but it does not endow any complexity that 
we have not seen before, 



D* 



St f Wt 



J^D^STdv 



1 



1 

of' 



-D* 



Wt 



1 



1 r Wt 



aT aT 



Wt 



1 



If we bring together the previous partial results we will obtain the expression, 

.-rT 



E 



SlaT \aT ^ a 



^{St) 



If we compare it with (0), we find the following relationship between V and F: 

V 



SlaT 



(11) 



(12) 



Since we have indeed closed expressions for all the Greeks, we may easily check the correctness 
of the above statements. We shall recover property ( |l^ ) of the European-style options in the 
next section. The above identities are very well known by practitioners although their proofs do 
not usually recall the integration by parts formula in the form we have introduced it here. 



4.2 The explicit computation 

The reason for the existence of such expressions for the Greeks of European-type options is due 
to the fact that there is a closed and tractable formula for the probability density function of 
St- This is the lognormal distribution that is written as 

p{x) = exp{-[log(x/go) - i,Tf/2a''T}. 
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When p{x) is available we can face the problem from a different perspective. In this case we are 
able to compute all the partial derivatives, starting from the explicit formulation for the price 
of the option, V, 



/•oo 

Jo 



x)p{x)dx, 



(13) 



usually just a formal expression, which now becomes handy. We can show this computing the 
value of A, in terms of partial derivatives of p{x): 



A 



— / e-'''^<^{x)p{x)dx = 



dp{x) 



[•oo 

Jo 



dSi 







d\ogp{x) 
dSo 



p{x)dx. 



Note that we get an expression that can be easily rewritten in a way that apparently resembles 
our previous results, since one integration by parts has been implicitly done, and a kernel 
naturally appears, 



A = S 



rT 



T 



dlogp{x) 



(14) 



'0 / a;=5j,J 

But we have not yet exploited the information we have about the functional form of p(x), 
which leads us to the same expression we have already obtain by means of Malliavin Calculus: 



E 



SqgT 

A similar procedure applies to the other Greeks. We will obtain Vega just replacing the 5*0 



with a (T in equation (14), 



V = E 



d\ogp{x) 



x=Sj 



and, after straightforward computations, we recover equation (10). The case of Gamma leads 
to an expression with a very similar flavor to what we have already seen. 



E 



-rT 



$(5' 



T 



d\ogp{x] 
dSo 



+ 



9^ logp{x) 



x=St 



that yields, again, the same result presented in (|lT|). We find therefore in this frame that the 
property stated in ([T^ ) is fulfilled by Vega and Gamma. 

We can then conclude that when we deal with European-style options, the Malliavin-related 
procedures presented above are equivalent to the result we attain if we directly differenciate the 
probability density function. 
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4.3 The vanilla options 



Besides the formal comparison with the previous case, the fact of knowing p(x) allows us, in 
principle, to completely compute all the Greeks once a payoff function has been selected. One 
of the most popular choice is the European, or vanilla, call whose payoff reads, 



$(X) = (X-K)+. 

Then can be easily derived the following expressions for the Greeks we have presented: 

1 rdi(K) 



(15) 



A 



V 

r 



V 27r 

1 ^-\di(K)]y2. 



where 



di{x) 



S'oV27r(j2r 



log(5o/x) + (r + -a^)T 



as it can be found in any textbook on financial derivatives Q. In conclussion, we are able 
to compute the different Greeks using the Malliavin-related formulas, and compare them with 
their theoretical values. We present in Fig. |l] and |2| the result of this procedure, for a given 
set of parameters, after Monte Carlo simulation. Only A and V are shown, since T would just 
be a replica of the second, due to equation (|T^). These examples show us how the outcome 
of the simulation progressively attains their own theoretical value, whereas the statistical error 
reduces. We notice however that the use of what we have labeled as "direct method", just 
performing Monte Carlo simulations starting from the rhs expression in (Q), would lead to an 
estimator with smaller variance, and therefore a better estimation. Those estimations do not 
appear in the Figures, in the sake of clarity. But we must remember that this technique can only 
be applied when the payoff is smooth enough. In our case, when payoff follows (^), Gamma 
cannot be computed in this way. 



5 The Asian-style options 

Here one considers the Greeks for options written on the average of the stock price y Sgds, 
instead of the final value St, as in European options. Note that in this particular case the 
density function of the random variable does not have a known closed formula. Delta in this 
case is given by 



_d_ 

dSi 



■E 







-rTf 



^ ( 4 / S,ds 







-rT 



So 



-E 



1 



^'[^ SAs ^ 







1 



Sudu 







There are various ways of doing the integration by parts. In the already cited literature we find 
in JtI] the following expression: 



-rT 



So 



-E 



1 



T 



^ I ^ / Ssds 



'2j^StdWt 
. <y !o Stdt 



+ 1 
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CO 
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Figure 1: Estimated value of Delta for an European call with parameters r = 0.1, a = 0.2, 
T = 1.0 (in years) and Sq = K = 100 (in arbitrary cash units), using Monte Carlo techniques. 



whereas a close variant of it, which involves (IBf), can be found in 



2e- 



-rT 



-E 



$1^^ Ssds 



St — Sq 
lo Stdt 



Of course, we may also use the same approach we have present in the previous sections, and 
obtain a third one: 



-rT 



A 



5*0 



-E 



where 



<S> 



1 (Wt <S^> 

< S > \ (T < S > 



1 



<s^> 



tStdt 
Jq Sydv 
llfStdt 



and 



are something similar to a first two moments. 

Although the two first expressions for A are statistically identical, their particular realiza- 
tions when perfoming numerical computation will slightly differ, even though the same series 
of random numbers is used. The last formula is definitely a brand new estimator with its own 
properties, among them its smaller variance is perhaps the most relevant one. We can observe 
these features in Fig. |3|, where we show the outcome of the Monte Carlo simulation using the 



three alternative instances. We have chosen again the functional form in (15) for the payoff, and 
the rest of parameters takes the same value we used in the making of the previous plots. 
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Figure 2: Estimated value of Vega for an European call with parameters r = 0.1, a = 0.2, 
T = 1.0 (in years) and Sq = K = 100 (in arbitrary cash units), using Monte Carlo techniques. 



Then, not all these formulas coincide and in fact, contrary to what is claimed in ^ there 
is no way to obtain the integration by parts that provides the minimal variance. The main 
reason being that this is equivalent to know the probability density of the random variable in 
question. To expose the main ideas that also appear in |^ one can note first that there is an 
integration by parts that is the "most" straightforward but highly unrealistic. For this, consider 
the generalized problem 



E 



S^ds 



T 



S^ds 



^'{x)xp{x)dx. 



Here p denotes the density of J^^ Sgds which exists and is smooth (it is an interesting exercise of 
Malliavin Calculus). Therefore one can perform the integration by parts directly in the above 
formula thus obtaining that 



E 



S^ds 



S^ds 



^{x){p{x) + xp'{x))dx 



E 



T 



S^ds 



1 + 



Jq Ssdsp'iJ^ Ssds) 



Now we procceed to prove that the above gives the minimal integration by parts in the sense of 
variance. Obviously it is not possible to carry out the simulations unless p' and p are known. 
Let us construct the set of all possible integration by parts. Suppose that y is a random variable 
such that 



E 



T 



S^ds 



S^ds 



E 



S,ds Y 
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Figure 3: Computed value of Delta for an Asian call with parameters r = 0.1, a = 0.2, T = 1.0 
(in years) and Sq = K = 100 (in arbitrary cash units), using Monte Carlo techniques, for the 
estimators presented in the main text. We have broken the interval of integration in 252 pieces, 
representing the approximate number of trading days in a year. The exact result turns out to 
be near the bottom of the graph, at about 0.65. 



for any function $ G Cp°°, then it is not difficult to deduce that 









E 













Here cr{x) denotes the u-algebra generated by x, and is the conditional expectation. 

Therefore the set of all possible integration by parts can be characterized as 



inf E 



T \ ^ 

S,ds] y2 











E 




= E 

















S.^ds] ■ (y 



7W = |y G L2(17); e 
Next in order we want to find the element in Y that minimizes 

As in [HI is not difficult to see which Y achieves the minimum. This is done as follows: 

^ _ Jq Ssdsp'iJ^ S,ds) \ 
pHo Ssds) J 







+E 



J \ p(fls.,ds) 
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since the mixed product is 0, due to the property of the set Ai. Therefore the minimum is 



is unknown in the case of Asian options. Therefore it is still an open problem to devise good 
ways to perform an efficient integration by parts so that the variance is made small rapidly and 
efficiently. 
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achieved 




This is clearly impossible to write explicitely as p 
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